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The complete 5.4 kton MINOS far detector has been taking data since the beginning of August 
2003 at a depth of 2070 meters water-equivalent in the Soudan mine, Minnesota. This paper presents 
the first MINOS observations of and charged-current atmospheric neutrino interactions based 
on an exposure of 418 days. The ratio of upward to downward-going events in the data is compared 
to the Monte Carlo expectation in the absence of neutrino oscillations giving: 

i?ur/down/<p/down = 0.62^“.I® (stat.) ± 0.02(sys.). 

An extended maximum likelihood analysis of the observed LjE distributions excludes the null 
hypothesis of no neutrino oscillations at the 98 % confidence level. Using the curvature of the 
observed muons in the 1.3 T MINOS magnetic field r-p and interactions are separated. The ratio 
of Fp to r'p events in the data is compared to the Monte Carlo expectation assuming neutrinos and 
anti-neutrinos oscillate in same manner giving: 

4;%p/-RFpVp = 0.96t“;“(stat.)±0.15(sj/s.), 

where the errors are the statistical and systematic uncertainties. Although the statistics are limited, 
this is the first direct observation of atmospheric neutrino interactions separately for and 77^. 
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I. INTRODUCTION 

Over the course of the past ten years the deficit 
of muon neutrinos from cosmic-ray showers in the at¬ 
mosphere has been firnfiy established by the Super- 
Kamiokande experiment 0 0) S 1^ and confirmed 
by the MACRO 0 and Soudan 201 experiments. The 
favoured interpretation of the data is neutrino 

oscillations. Recent results from the Super-Kamiokande 
experimental provide direct evidence for atmospheric 
neutrino oscillations and yield best fit oscillation pa¬ 
rameters of (Am 23 , sin^ 26 * 23 ) = (0.0024eV^, 1 . 0 ), where 

Am^g = Iwg—TO 2 I. Results from the K2K experiment ^3 

provide further confirmation of the oscillation 

hypothesis. 

The 5.4kiloton (kt) mass of the recently constructed 
MINOS (Main Injector Neutrino Oscillation Search) far 
detector is much less than the ~ 25 kt fiducial mass of 
the Super-Kamiokande detector. However, it does pos¬ 
sess one unique advantage, namely it is the first large 
deep underground detector to have a magnetic field. This 
allows studies of neutrino flavour oscillations for neutri¬ 
nos and anti-neutrinos separately. A separate measure¬ 
ment of and oscillations could provide constraints 
on CPT violating models [T^IT^ I which have been invoked 
to accommodate simultaneously the solar, atmospheric 
and LSND[i| neutrino oscillation data. It should be 
noted that a number of recent studies have indicated 
difficulties with the CPT violating models (see for ex¬ 
ample M)- Nevertheless, a direct measurement of 
and oscillations is of interest. In addition, MINOS is 
unique in its ability to provide an accurate measurement 
of the neutrino energy and direction for all contained- 
vertex charged-current (CC) interactions. 

This paper presents first results on atmospheric neutri¬ 
nos from the MINOS experiment. Here, only results from 
CC events with neutrino interaction vertices con¬ 
tained inside the detector volume are considered; results 
from events where the neutrino interacts in the surround¬ 
ing rock will be the subject of a separate publication. The 
data used were recorded between August 2003 and Febru¬ 
ary 2005 and correspond to a livetime of 418 days giving 
an exposure of 6.18 kiloton-years (4.54 kiloton-years fidu¬ 
cial). The data are compared to the expectation in the 
absence of neutrino oscillations and the favoured hypoth¬ 
esis of <-> Vr oscillations with Am^g = 0.0024eV^ and 
sin^ 2023 = 1-0. The first direct results showing charge 
separated and 77^ atmospheric neutrino interactions 
are presented. 

II. THE MINOS DETECTOR 

The MINOS far detector is located at a depth of 
2070 meters-water-equivalent (mwe) in the Soudan mine. 
Northern Minnesota. The far detector is a steel- 
scintillator sampling calorimeter consisting of two super¬ 
modules (SM) separated by a gap of 1.1m. The detector 


consists of octagonal planes of 2.54 cm thick steel followed 
by planes of 1 cm thick extruded polystyrene scintillator 
and a 2 cm wide air gap. The first and second SMs are 
comprised of 248 and 236 scintillator planes respectively. 
Each SM is magnetized to an average value of 1.3 T by a 
15 kA current loop which runs through the coil hole along 
the detector central axis and returns below the detector. 
Each scintillator plane is made up of 192 strips of width 
4.1cm and of length between 3.4 — 8.0 m depending on 
position in the plane. The strips in alternating planes 
are oriented at ±45° to the vertical thereby providing 
two orthogonal coordinates^. The scintillation light is 
collected using wavelength shifting (WLS) fibers embed¬ 
ded within the scintillator strips. The WLS fibers are 
coupled to clear optical fibers at both ends of a strip 
and are read out using 16-pixel multi-anode photomulti¬ 
plier tubes (PMTs). The signals from eight strips, sep¬ 
arated by approximately 1 m within the same plane, are 
optically summed (multiplexed) and read out by a single 
PMT pixel. The multiplexing pattern is different for the 
two sides of the detector, which, for a single hit, enables 
the resulting eightfold ambiguity to be resolved. For all 
types of event, the ambiguities are efficiently resolved in 
software using additional information from timing and 
event topology. 

The detector is optimized for detecting beam neutri¬ 
nos coming from the direction of Fermilab. For the study 
of atmospheric neutrinos the planar structure presents a 
particular problem: cosmic-ray muons traveling almost 
parallel to the scintillator planes can penetrate deep into 
the detector by traveling in the steel or air between the 
planes. To reject this source of background a scintillator 
veto shield surrounds the upper part of the main detec¬ 
tor. The veto shield is constructed from the same scintil¬ 
lator modules as used in the main detector but with the 
orientation of strips aligned along the z-axis. The veto 
shield comprises a “ceiling” section above the detector, 
consisting of two scintillator layers, and “wall” sections 
along each of the two sides of the detector formed from 
a single scintillator layer. 


A. Data Acquisition and Trigger 

The output signals from each PMT pixel are digitized 
and time-stamped (with a 1.5625 ns precision) by the 
VME-based front-end electronics. The signals from the 
pixels are digitized by 14-bit analogue-to-digital convert¬ 
ers (ADC) when the dynode signal from the PMT ex¬ 
ceeds a programmable threshold, corresponding to ap- 


^ The MINOS right-handed coordinate system has the 2 :-axis de¬ 
fined along the detector axis pointing away from Fermilab and 
the j/-axis vertical. The alternating scintillator planes provide 
measurements of the U and V coordinates which are related to 
X and y hy U = ^{x + y) and V = -^{y - x). 
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proximately one third of a photo-electron. To reduce the 
data flow, the pedestal corrected signals are only writ¬ 
ten to the DAQ output buffers if two out of thirty-six 
channels on the same readout board are above thresh¬ 
old. These thirty-six channels correspond to the readout 
on one side of the detector from a contiguous group of 
either 20 or 24 planes. The raw data rate is approxi¬ 
mately 8MBs“^. The raw data are transferred to a PC 
based trigger farm where the data are divided into blocks 
bounded by regions of 100 clock ticks (156 ns) or more 
where no detector activity has been recorded. The pri¬ 
mary trigger algorithm, applied to these blocks of data, 
requires there to be activity in at least 4 planes out of 
any contiguous group of 5 planes. The veto shield is read 
out in the same manner as the main detector except that 
the two out of thirty-six requirement is not applied and 
the dynode threshold is set to a level corresponding to 
approximately one and a half photo-electrons. 

The MINOS far detector front-end electronics and data 
acquisition system are described in detail in 0 and da 

B. Detector Calibration 

A minimum ionizing particle crossing at normal inci¬ 
dence to a plane gives a combined signal of approximately 
10 photo-electrons (PEs) registered by the PMTs at the 
two ends of the strip. The detector is calibrated using 
both a dedicated LED system0| and cosmic-ray muons. 
The ADC to PE calibration is performed using the LED 
system and the cosmic-ray muon sample is then used to 
give a uniform response across the detector. From stud¬ 
ies of cosmic-ray muons in the MINOS detector 0|, the 
current uncertainty in the PE calibration is 5 %. Cosmic- 
ray muons are also used to calibrate the recorded times. 
After calibration, a single hit timing resolution of approx¬ 
imately 2.3 ns is achieved. The timing calibration tracks 
all hardware changes. 

III. DATA AND MONTE CARLO 

The data described in this paper were recorded in the 
18 month period from August 2003 to February 2005. 
Only data taken when the MINOS far detector, including 
the veto shield, was fully operational are used. The final 
data sample corresponds to a livetime of 418 days giv¬ 
ing an exposure of 6.18 kiloton-years (4.54 kiloton-years 
fiducial). 

The selection of contained-vertex neutrino interactions 
was optimized using a GEANT3|0 simulation of the 
MINOS detector. For the simulation of atmospheric neu¬ 
trino events the 3D flux calculation of Barr et al. was 
used (Bartol 3D). The NEUGEN3 programwas used 
to simulate the neutrino interactions (cross sections and 
hadronic final states). The earlier ID flux calculation 
from the Bartol groupj^ (Bartol ID) and the 3D cal¬ 
culation of Battistoni et were used to assign sys¬ 


tematic uncertainties. The response of the MINOS de¬ 
tector to electrons, muons and hadrons has been studied 
in a test beam at the CERN PS using the 12.5 ton MI¬ 
NOS calibration detector[0. The test beam detector 
was constructed and read out in the same manner as 
the MINOS far detector. The interactions of hadronic 
particles are modeled with the GCALOR package 0. 
which is found to give a reasonable description of low 
energy hadronic interactions in the MINOS calibration 
detector|^, rather than the default version of GHEISHA 
(see and references therein). The “SLAG version” 
of GHEISHA, which also provides a reasonable descrip¬ 
tion of the test beam data, is used as an alternative model 
for hadronic interactions. A Monte Garlo (MC) sample of 
atmospheric neutrino interactions corresponding to over 
1000 live-years was generated and used to optimize both 
the reconstruction algorithms and the event selection cri¬ 
teria. Two large cosmic-ray muon background samples 
were generated: a sample of 19 million events full spec¬ 
trum (corresponding to approximately 280 days livetime) 
and a further 2 million events with E^ < 2 GeV (cor¬ 
responding to a livetime of 4.1 years) as lower energy 
events are an important component of the cosmic-ray 
muon background to the contained-vertex atmospheric 
neutrino selection. A 10 % uncertainty on the normalisa¬ 
tion of the cosmic-ray background is assigned. The error 
reflects the different normalisation obtained when nor¬ 
malising to the entire cosmic-ray sample or normalising 
to just those cosmic-ray muons which stop in the detector 
(these form the main background to the event selection 
described below). It should be noted that for the results 
presented in this paper the cosmic-ray background in the 
selected event sample is estimated from data. The 10 % 
uncertainty in the cosmic-ray normalisation is only used 
when comparing data and Monte Garlo at various stages 
in the event selection. 


A. Flux Normalization and Systematic 
Uncertainties 

The theoretical prediction for the atmospheric neu¬ 
trino event rate has large uncertainties from the primary 
cosmic-ray flux, hadron production models and neutrino 
interaction cross sections. The analysis of the Soudan 2 
dataQ indicates that the combined pred iction of 
the Bartol 3D modelj^ and the NEIJGEN3|22l| neutrino 
cross section model should be scaled by 0.88 ± 0.070, 
where the error is statistical and it is implicitly as¬ 
sumed that atmospheric electron neutrinos are not oscil¬ 
lating. This normalization result is compatible with the 
results from a fit to the Soudan 2 data including oscilla¬ 
tions. The MINOS and Soudan 2 detectors are located 
in the same mine (i.e. at the same geomagnetic lati¬ 
tude) and both are constructed from steel. Consequently, 
for the analysis presented here the Soudan 2 scale factor 
0.88 ± 0.07 is used to correct the combined event rate 
predictions from the Bartol 3D flux model and the NEU- 
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GENS neutrino interaction model. An additional 5 % un¬ 
certainty is added in quadrature to that estimated by the 
Soudan 2 collaboration to account for differences arising 
from the different energy thresholds (300 MeV in the case 
of Soudan 2 compared to ^ 500 MeV for MINOS). Finally 
an additional 2.5 % uncertainty is assigned to account for 
the different phases in the solar cycle for the Soudan 2 
and MINOS data sets[^; because of the relatively high 
neutrino energy threshold the selected atmospheric neu¬ 
trino rate in MINOS depends only weakly on the phase 
in the solar cycle. The resulting total systematic uncer¬ 
tainty on the expected event rate is estimated to be 10 %. 


IV. EVENT RECONSTRUCTION 

The MINOS detector is optimized for beam neutrinos 
originating from Fermilab. Due to the curvature of the 
Earth, beam neutrinos enter the detector from below the 
horizontal at an angle of 3.3° with respect to the 2 ;-axis. 
The standard MINOS reconstruction software has been 
developed for these events. The analysis presented here 
uses reconstruction software optimized for atmospheric 
neutrinos!^. 

The first stage of the event reconstruction removes 
the eightfold ambiguity in the association of raw hits 
to strips. This is performed utilizing information from 
both strip ends. For cosmic-ray muons, an average of 
99 % of the recorded pulse height is associated with the 
correct strip. At this stage the data are in the form 
of two 2D event views U — z and V — z. An example 
event display of a cosmic-ray muon is shown in Figure ^ 
Tracks and showers are reconstructed independently in 
each view; the two views are then matched to obtain 
a three-dimensional event. For cosmic-ray events that 
leave hits in both the veto shield and main MINOS far 
detector, the root-mean-square (rms) difference in times 
recorded in veto shield and the detector is 4 ns, allowing 
association of veto shield hits (indicated in Figure to 
activity in the main detector. 

A charged-current muon neutrino event is in general 
reconstructed as a muon track and a hadronic shower. A 
typical 1 GeV muon will traverse approximately 25 planes 
at normal incidence. Reconstructed tracks are required 
to consist of at least 8 planes (corresponding to a min¬ 
imum energy of 0.4GeV). For muons which start and 
stop within the detector volume the muon momentum 
is determined from range with a resolution of approxi¬ 
mately (cTp/p)^ = 0.06^ -I- (0.045/p)^ for muons traveling 
at normal incidence to the detector planes (where p is 
measured in GeV/c). The first term is dominated by 
fluctuations in energy loss and the second is dominated 
by sampling. For events where the muon exits the detec¬ 
tor, the muon momentum is obtained from the curvature 
of the track in the magnetic field. For the selected GG 
atmospheric interactions, where the momentum is 

determined from curvature, the average momentum res¬ 
olution is approximately = [0.1^ -I- (0.3/p)^] GeV“^ 




X / m 


FIG. 1: An example of a cosmic-ray muon event in the MI¬ 
NOS far detector. The detector readout corresponds to the 
two orthogonal U — z and V — z views. The size of the points 
gives an indication of the pulse height for each scintillator 
strip hit. The large dark points shown in the x — y view 
indicate in-time activity in the veto shield. 

(where p is measured in GeV/c). However, the resolu¬ 
tion obtained from individual events depends strongly 
on how much of the trajectory of the muon is observed 
before it exits the detector and on the orientation of 
the trajectory relative to the local magnetic field. The 
hadronic energy is obtained by summing the pulse height 
in a shower which is spatially associated with the start 
of the track. The energy scale is obtained from Monte 
Garlo using the GGALOR model of hadronic show¬ 
ers, which from the test beam results is found to provide 
a good description of the detector response to single 
and protons 1^. The hadronic energy resolution is ap¬ 
proximately oe/E ^ 0.55/-\/F, where E is measured in 
GeV. 

For the study of atmospheric neutrinos it is necessary 
to determine whether the reconstructed track is upward 
or downward going. A relativistic normal incidence par¬ 
ticle traverses ten planes in approximately 2 ns which, 
when compared to the single hit resolution of 2.3 ns, is 
sufficient to identify the direction of most selected events 
with little ambiguity. The sense of the direction of muon 
tracks is determined by comparing the hit times along the 
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FIG. 2: The efficiency for correctly reconstructing stopping 
muon events as downward-going as a function of number of 
planes in the reconstructed track. 


reconstructed track with the hypotheses that it is either 
upward or downward going (assuming that the particle 
is traveling at the speed of light). The rms deviations 
of hit times about each of the two hypotheses are cal¬ 
culated, RMSup and RMSdown- The hypothesis with 
the smallest rms is chosen. In addition, the magnitude 
of RMSup — RMSdown provides a measure of the qual¬ 
ity of the direction determination. To test the perfor¬ 
mance of the algorithm a sample of stopping cosmic-ray 
muons is used (all of which are traveling downward). Fig¬ 
ure [3 shows a comparison of the data and Monte Carlo 
efficiencies for correctly identifying a stopping muon as 
downward-going as a function of the number of planes 
the track crosses. The average efficiency is above 94% 
for even the shortest tracks and rapidly increases to bet¬ 
ter than 99% for events with hits in 12 planes. The 
efficiency in data agrees with that from Monte Carlo to 
better than 1 %. 

The curvature of /fj,~ tracks in the magnetic field al¬ 
lows the charge sign to be determined. FigureElshows the 
distribution of the reconstructed charge divided by mo¬ 
mentum, Q/p, divided by its error, for cosmic-ray muons 
that stop in the detector. Two peaks, corresponding to 
fj,~ and events, are clearly seen. The widths of the two 
peaks in data and MC agree to better than 2.5%. For 
the event samples considered here, the charge is 

cleanly identified over the approximate momentum range 
1 — 10 GeV. The efficiency decreases for low momentum 
tracks due to the limited number of planes crossed. For 
high momentum tracks, which typically leave the detec¬ 
tor, the charge identification efficiency decreases as only 
the limited curvature at the start of the track is mea¬ 
sured. 



FIG. 3: The reconstructed distribution of (Q/p)/cr(Q/p) for 
stopping muon events in data and Monte Carlo. 


V. EVENT SELECTION 

At a depth of 2070 mwe the cosmic-ray muon rate is 
approximately 50000 events per day in the MINOS de¬ 
tector. This rate should be compared to the expected 
signal rate of 0.54 ± 0.05 atmospheric CC interac¬ 

tions per day^, where the uncertainty is from the 10% 
uncertainty in the expected event rate (discussed in Sec¬ 
tion IIII All . In order to achieve a signal-to-background 
ratio of ten-to-one it is necessary to identify the signal 
events efficiently whilst reducing the background by a 
factor of 10®. The event selection is designed to identify 
both fully-contained (FC) and partially-contained (PC) 
events. In FC events the entire event is contained 
within the fiducial volume. In PC events the neutrino 
vertex is within the fiducial volume but the produced 
muon exits the detector. 


A. Preselection 

Candidate CC neutrino interactions are required 
to have a reconstructed track passing some basic quality 
requirements. The majority of the background is rejected 
by event containment requirements which are applied at 
both the hit and reconstructed track level. The sense of 
the track direction (up/down) is determined from timing 
as described previously. The start of the track, which 
is considered to be the neutrino interaction vertex, is 
required to lie within the detector fiducial volume. The 
fiducial volume is defined as the octagonal region which 
is at least 50 cm from the detector edges in the xy plane 
and at least 5 planes from the start and end of either SM. 


^ The signal rate of 0.54 ± 0.05 (no oscillations) corresponds to 
I'll I I'll CC interactions where the muon deposits energy in at 
least eight planes (before fiducial cuts). 
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In addition, the region within 40 cm of the axis of the coil 
hole, which has a diameter of 25 cm, is excluded from the 
fiducial volume. This cut is enlarged to 1 m in the first 
and last 10 planes of the detector. The event sample is 
sub-divided into FC and PC events depending on whether 
the end of the track also lies within the fiducial region. 

Event containment cuts are also made at the hit level 
to reduce the sensitivity to possible reconstruction errors 
where not all hits are correctly associated to the recon¬ 
structed track. For this purpose, the fiducial volume re¬ 
quirement of 50 cm from the detector edges is relaxed to 
30 cm. In order to apply the containment cuts at the hit 
level it is necessary to convert the two-dimensional coor¬ 
dinates of a single hit into a point in space. This conver¬ 
sion is achieved by using the mean value of the orthog¬ 
onal (U/V) coordinate in the surrounding two planes. 
Hits outside the fiducial volume are then assigned to 
the nearest octagonal edge/SM end. Edges/ends with 
summed pulse height equivalent to more than six PEs 
outside the fiducial region are tagged as being uncon¬ 
tained. Candidate FC (PC) events are required to have 
no (one) such region. The containment cuts reject ap¬ 
proximately 99.9% of the cosmic-ray background whilst 
retaining 77% of CC Vfilv^ interactions in the detector 
volume. The inefficiency for signal events is primarily a 
fiducial effect; the containment cuts retain 99 % of CC 
interactions in the fiducial region which produce a 
muon which spans at least six planes. 

Candidate FC events are required to have a recon¬ 
structed track consisting of hits in at least eight planes. 
The PC event selection criteria are optimized separately 
for upward and downward-going events as the back¬ 
grounds for the two categories are very different. To 
ensure the track direction is well determined, candidate 
PC events are required to have a track of at least 1 m in 
length and which consists of hits in at least ten planes. 

B. Fully Contained and Downward Partially 
Contained Event Selection 

The dominant backgrounds in the FC and downward 
PC samples arise from steep cosmic-ray muons which en¬ 
ter the detector at small angles to the detector planes. 
By traveling in the steel or air between the scintillator 
planes, such events can penetrate a significant distance 
into the fiducial volume before leaving a detectable sig¬ 
nal. The selection of FC and downward-traveling PC 
CC iZ/i/P/j interactions aims to greatly reduce this back¬ 
ground and proceeds in four stages: 

i) Cosmic-Ray Rejection (Trace cut): The re¬ 
constructed track is extrapolated back to the out¬ 
side of the detector and the distance traversed in 
the direction perpendicular to the detector planes 
is calculated, A^. Events with small values of 
correspond to steep tracks which when extrapo¬ 
lated to the detector edge traverse only a few scin¬ 
tillator planes. Figure^lshows the A^ distribution 



FIG. 4: The reconstructed Az distribution for events pass¬ 
ing the containment cuts. The solid histogram indicates the 
MC expectation for cosmic-ray background. The points with 
error bars show the observed data. The dashed histogram 
indicates the expected distribution for atmospheric neutrino 
events (without oscillations). 

for MC cosmic-ray muons and CC interac¬ 

tions. Events are rejected if A^ < 0.5 m. Eigure^l 
also shows the A^ distribution for data which is in 
reasonable agreement with the MC expectation. 

ii) Event Topology: About half of the remaining 
background consists of cosmic-ray muon tracks that 
bend in the magnetic field and turn over in the z di¬ 
rection. Such events will leave hits in two separate 
positions in a particular plane. In addition, these 
events typically have large pulse height in the plane 
where the muon turns around in z. This category 
of background event is rejected using the pulse- 
height weighted deviations of the hits in the U — z 
and V — z views from the fitted track. The pulse- 
height weighted mean, (Ayv)) and the pulse-height 
weighted rms deviation of hits from the track, 
(Auv) ^) are calculated. Events are rejected if there 
is large scatter about the track, (Auy )2 > 0.5 m, 
or if the pulse-height weighted mean deviation from 
the track lies significantly above the reconstructed 
track, (Auv) > 0.25 m. These empirically deter¬ 
mined cuts are applied separately to the hits in 
both the U — z and V — z views. In addition, the 
event vertex is defined as the first hit on the track 
taking the highest end (largest y) as the start of 
the track. The maximum displacement from the 
event vertex of the hit strips which lie within ±4 
planes of the event vertex is found, A^““. Events 
are rejected if A^““ > 1.25 m. 

hi) Vertex Pulse Height/Direction: After the 
topology cut, the signal-to-background ratio is ap¬ 
proximately 1:5. The remaining background con¬ 
sists of steep cosmic-ray muons which travel nearly 
parallel to the scintillator planes and therefore tend 






ad/‘0 ad/ 


to give a large pulse height signal in a single plane 
near the beginning of the track. These events are 
often poorly reconstructed due to the difficulties 
of reconstructing tracks for events at small an¬ 
gles to the detector planes. Figure El shows, for 
signal and background, the total pulse height in 
the event vertex region, Qvtx, plotted against the 
cosine of reconstructed zenith angle at the high¬ 
est end of the track^ and the modulus of recon¬ 
structed track direction cosine with respect to the 
z axis, I cos 9^ \. The vertex pulse height is defined 
as the maximum number of PEs observed in a single 
plane within ±4 planes of the event vertex (defined 
above). The background is characterized by being 
steep and having large Qvtx- Events are rejected if 
they have Qvtx > 300 PEs. Steep events, defined as 
having |cos0zen| > 0.7 and |cos0z| < 0.5, are re¬ 
quired to satisfy Qvtx < 100 PEs. The above event 
charge/direction cuts are not applied to events with 
track lengths of greater than 20 detector planes, as 
the steep background events tend to cross relatively 
few planes. 




® The zenith angle, 6zen, is defined as tt minus the angle between 
the reconstructed track direction and the local vertical (the y- 
axis). Negative values of the cosine of the zenith angle correspond 
to tracks which are reconstructed as upward-going. 


FIG. 5: The MC distributions of the vertex pulse height, 
Qvtx, plotted against the modulus of the cosine of the angle 
with respect to the detector z-axis, | cos 9z \ and the cosine 
of the reconstructed zenith angle, cos^zen. Plots are shown 
for both the cosmic-ray muon background and for the atmo¬ 
spheric neutrino signal for all events passing the containment 
cuts. The hatched areas represent the regions rejected by the 
“Vertex Pulse Height/Direction” cuts. 


iv) Veto Shield: The cuts listed above result in a 
signal-to-background ratio of approximately 1 : 
2. Additional background is removed by reject¬ 
ing events with activity in the veto shield within 
a ±100 ns window around the event time, resulting 
in a signal-to-background ratio of approximately 
20 : 1 . 


C. Upward Partially Contained Events 

The background to the upward-going PC event se¬ 
lection is dominated by cosmic-ray muons which stop 
in the detector and are reconstructed as upward rather 
than downward-going. The cuts to remove this source of 
background are based on timing information and iden¬ 
tify events which are unambiguously upward-going. The 
event selection cuts fall into two categories: 

i) Event Topology: For the upward PC selection, 

the relatively small number of badly reconstructed 
events passing the preselection are removed us¬ 
ing a subset of the topology cuts employed in the 
FC/downward PC analysis. Events are rejected if 
Qvtx > 300 PEs or > 1.25 m. 

ii) Track timing rms: The expected times of hits 
on the track are calculated for the hypotheses of 
an upward-going and a downward-going track. The 
rms scatter of the difference between observed and 
expected hit times for these two hypotheses is used 
to identify upward-going events. The event is re¬ 
quired to be consistent with the upward hypothe¬ 
sis and significantly more compatible with the up¬ 
ward hypothesis than the downward hypothesis: 
RMSup < 4.33 ns and RMSup ~ RMSdown < 
— 1.66 ns (these numbers should be compared to 
the single hit timing resolution of 2.3ns). FigureEl 
shows the effect of the main timing cut, namely 
RMSup — RMSdown, for data compared to the 
Monte Carlo expectation. 


D. Performance 

The event selection reduces the background from 
cosmic-ray muons by a factor of 4 x 10®. In Monte Carlo 
the efficiency for CC neutrino interactions where 

the interaction occurs within the fiducial region and the 
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RMSup-RMSpo^y^/ns 


FIG. 6: The distribution of RMSup — RMSdown for events 
passing all other cuts in the upward partially contained event 
selection. The data are shown by the points with error bars; 
the total Monte Carlo expectation is shown by the solid his¬ 
togram with the expected atmospheric neutrino contribution 
(no oscillation) shown by the dashed histogram. The cut is 
indicated by the arrow. 



FIG. 7: The expected number of CC events that inter¬ 

act within the detector fiducial volume at the different stages 
in the selection (418 days exposure). In addition to the ex¬ 
pected energy distribution before selection and the final se¬ 
lected distribution, the expected numbers of events are shown 
at several stages in the analysis: events which pass the trigger 
requirements; events for which there is a reconstructed track; 
and events which pass the preselection cuts. 


muon traverses eight or more scintillator planes is 70%. 
Figure [7| shows the expected energy distribution for CC 
for the various stages in the event selection. The 
effective lower limit on the neutrino energy of the selected 
events is approximately 0.5 GeV. For low energy CC neu¬ 
trino interactions the efficiency is low because tracks are 
only reconstructed if they span at least eight detector 
planes. 

The numbers of events surviving at different stages in 
event selections are listed in Tabled Reasonable agree¬ 


ment between data and Monte Carlo is seen at each stage. 
The final veto shield requirement rejects 149 events in 
data, consistent within one standard deviation with the 
MC expectation of 170±24, where the uncertainty is from 
normalization and MC statistics. For the results in this 
paper, the cosmic-ray muon background in the combined 
FC and downward PC event sample is estimated from 
data rather than relying on MC. From the 149 events re¬ 
jected by the veto shield cuts the remaining background 
is estimated to be 4.4 ± 0A{stat.) ± 0.3(sj/s.) events us¬ 
ing the veto shield efficiency of 97.1 ± 0.2% (described 
below) and taking account of the expected number of 
neutrino events rejected by the veto shield requirements. 


E. Event Selection Systematic Uncertainties 

The systematic uncertainties on the event selection effi¬ 
ciency and cosmic-ray muon background have been stud¬ 
ied in detail. In each case the impact of the systematic 
effect on the MC expectation for the number of selected 
events is estimated. In addition, because the selection is 
not up-down symmetric, systematic uncertainties are cal¬ 
culated for the MC expectation for the ratio of upward- 
going to downward-going events. The total uncertainty 
on the selection efficiency for atmospheric neutrino events 
is estimated to be 3.3 %. The contribution to the system¬ 
atic uncertainty on the up-down ratio from experimental 
effects is estimated to be 3.1%. The contributions to 
these systematic errors are discussed in detail below. 

Veto Shield: The efficiency of the veto shield cut is 
determined directly from data in two independent ways. 
Firstly, a sample of cosmic-ray muons that stop in the de¬ 
tector and have | cos0zen| > 0.5 is selected. These events 
occupy a similar region of phase space to the background. 
The veto shield cut rejects 97.06 ± 0.03% of this sam¬ 
ple. A second estimate of the veto shield efficiency is 
obtained by relaxing the event selection cuts until the 
sample is dominated by background {i.e. an expected 
signal fraction of less than 2 %). The veto shield cut 
rejects 96.2 ± 0.2% of this sample which, when the ex¬ 
pected signal (assuming Am^s = 0.0024eV^) is taken 
into account, leads to an estimated veto shield efficiency 
of 97.3 ± 0.2%. From these two tests the veto shield 
efficiency is estimated tobe97.1±0.2%, where the cen¬ 
tral value is taken from the high statistics stopping muon 
sample and a systematic error of 0.2 % is added reflecting 
the difference between the two methods. 

The fraction of signal events rejected due to accidental 
coincidences with hits in the veto shield is estimated by 
overlaying veto shield hits obtained from special min¬ 
imum bias data taking runs onto Monte Carlo atmo¬ 
spheric neutrino events. The estimated fraction of sig¬ 
nal events rejected due to spurious veto shield hits is 
2.2±0.4 %, where the error represents systematic time de¬ 
pendent variations. In addition, from Monte Carlo stud¬ 
ies it is estimated that 0.3 ± 0.1 % of the selected signal 
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Cuts 

Data 

Monte Carlo Expectation 



total 

cosmic muon 

CC 

Fully Contained and Downwards Partially Contained 

Preselection 

41571 

38253 ± 3987 

38121± 3987 

125 ± 13 

Trace 

1525 

1513 ± 153 

1395 ± 153 

112 ± 12 

Topology 

560 

494 ± 48 

384 ± 47 

104 ± 11 

Vertex/Direction 

243 

277 ± 26 

170 ± 24 

102 ± 11 

Veto Shield 

94 

110 ± 11 

4.9 ±0.7 

100 ± 10 

Upward Partially Contained 

Preselection 

427 

408 ± 47 

384 ± 47 

24 ± 2 

Topology 

364 

359 ± 42 

336 ± 42 

22 ±2 

Timing 

13 

18 ±2 

< 0.36 (68% C.L.) 

17 ± 2 


TABLE I: The numbers of data events after the different stages of the event selection compared with the MC expectation 
from cosmic-ray background events and CC atmospheric and 77^ events. The atmospheric neutrino numbers are the MC 
expectations for no oscillations. The uncertainties include MC statistics and systematic uncertainty on the normalization 
(±10% for the cosmic-ray background and 10% for the atmospheric neutrino sample) and a 3.3% systematic uncertainty 
on the selection efficiency for CC atmospheric and 77^ events. The numbers in the total column include neutral current 
interactions, v^jVe. CC interactions and interactions of neutrinos in the surrounding rock. 


downward-going events will be rejected due to hits in the 
veto shield associated with the neutrino interaction. 

Hadronic Response: The event selection efficiency de¬ 
pends on the detector response to hadrons and conse¬ 
quently the hadronic interaction model used. Compar¬ 
isons of GCALOR and GHEISHA show no evidence for 
any significant difference in overall selection efficiency or 
reconstructed up-down ratio. The GCALOR model is 
found to provide a good description of the response of 
the detector to single hadrons. Systematic errors of 2.5 % 
on the selection efficiency and 3.0% on the up-down ra¬ 
tio are assigned; in both cases the estimates reflect the 
Monte Carlo statistical precision of the comparison. 

Scintillator Light Calibration: The overall calibra¬ 
tion of the MINOS far detector is currently known to 
5 %. The MC response is tuned to agree with cosmic- 
ray muon data and has a corresponding 5 % uncertainty. 
Because the selection cuts use pulse height information, 
this leads to systematic errors of 0.6 % on the selection 
efficiency and 0.3% on the up-down ratio. 

Timing Calibration/Resolution: The timing calibra¬ 
tion for each scintillator strip is determined from data 
in a manner that tracks hardware changes. The uncer¬ 
tainty on the timing calibration for the individual strips is 
0.3 ns, i.e. significantly less than the single hit resolution 
of 2.3 ns. The effect on the selection is negligible. A more 
significant effect is that the single hit resolution in MC is 
better than that in the data, 2.2 ns compared to 2.3 ns. 
This is due to an incomplete simulation of the electronics 
readout. For this reason the times of the Monte Carlo hits 
are smeared by a Gaussian of width 0.7 ns. The difference 
between the the selection efficiencies before and after this 
smearing are compared. The overall selection efficiency 
for the smeared MC is reduced by 1.0%. The effect on 
the up-down ratio is small (0.1 %). These differences are 
used as estimates of the systematic uncertainties. 


Muon d£!/dx: One of the main cuts in the event selec¬ 
tion is the requirement that tracks leave hits in at least 
8 scintillator planes. Consequently, the event selection 
efficiency is sensitive to the Monte Carlo simulation of 
muon energy loss. The simulation of muon energy loss 
depends on the underlying simulation of the physics pro¬ 
cesses and the knowledge of the chemical composition of 
the MINOS detector. An uncertainty of 3 % in the muon 
range is assumed. These uncertainties result in system¬ 
atic uncertainties of 1.7% on the selection efficiency and 
0.3 % on the up-down ratio. 

Neutron Background: The background from cosmic- 
ray induced neutrons has been studied using a GEANT 
4 simulation|32| of muon nuclear interactions in the rock 
and is found to be negligible. In a Monte Garlo sample 
corresponding to four times the data exposure no neutron 
event passed even the early stages of the event selection. 

VI. RESULTS 

From the 418 days exposure considered in this pa¬ 
per, 107 candidate contained events are selected. The 
107 selected events are consistent with both the expec¬ 
tation of 127 ± 13 events assuming no neutrino oscilla¬ 
tions, and with the expectation of 96 ± 10 events assum¬ 
ing Amlg = 0.0024eV^ and sin^ 2023 = 1-0. The back¬ 
ground contribution from cosmic-ray muons, 4.4 ± 0.5, is 
obtained from data as described above. In addition, there 
is an expected background of 4.5 ±0.5 from the combina¬ 
tion of neutral current interactions and Vejve GG interac¬ 
tions. The error in the MG expectation is dominated by 
the uncertainty on the neutrino flux x interaction cross 
section which is estimated to be 10 %. Table im gives a 
breakdown of the various contributions to the expected 
event rates. 

The xy positions of the reconstructed neutrino inter¬ 
action vertices is shown in Figure |H1 There is no evidence 
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Selection 

Data 

cosmic fi 

Expe 

CC 

ctation (no 
Ve/l'e CC 

oscillatio 

NC 

ns) 

Rock Vfj, 

j CC 

FC 

69 

3.9 ±0.4 

81.2 ±8.5 

2.5 ±0.3 

2.0 ±0.2 

0.3 ±0.1 

- 

PC Down 

25 

0.6 ±0.2 

18.5 ± 1.9 

0.1 

— 

0.1 

— 

PC Up 

13 

< 0.36 

17.4 ± 1.8 

- 

- 

0.1 

- 

Total 

107 

4.4 ±0.5 

117.1 ± 12.2 

2.6 ±0.3 

2.0 ±0.2 

0.5 ±0.1 

- 




Expectation {Ami 

3 = 0.0024 eV“) 


FC 

69 

3.9 ±0.4 

58.4 ±6.1 

2.5 ±0.3 

2.0 ±0.2 

0.2 

0.7 ±0.1 

PC Down 

25 

0.6 ±0.2 

17.5 ± 1.8 

0.1 

— 

0.1 

— 

PC Up 

13 

< 0.36 

9.2 ± 1.0 

- 

- 

0.1 

0.5 ±0.1 

Total 

107 

4.4 ±0.5 

85.1 ±8.9 

2.6 ±0.3 

2.0 ±0.2 

0.4 ±0.1 

1.2 ±0.1 


TABLE II: The numbers of data events in each selection category compared to the expectation from different sources. The 
MC expectations from neutrino interactions are given for both no oscillations and sin^ 2023 = 1-0 and Am |3 = 0.0024eV^. 
The column referring to “rock refers to muons which are produced by neutrino interactions in the surrounding rock. The 
cosmic muon backgrounds in the FC and PC Down samples are estimated from data events passing all selection cuts with the 
exception of the veto shield. Entries marked as — indicate expectations of less than 0.05 of an event. For the entries where no 
error is quoted the error is less than 0.05 of an event. 



x/m 


FIG. 8: The reconstructed x — y positions of the neutrino 
interaction vertices for the 107 selected events. The vertex 
is defined as the start of the track, which is determined from 
timing. The solid lines indicate the active region of the MI¬ 
NOS detector and the dotted lines indicate the boundaries of 
the fiducial volume. 



FIG. 9: The reconstructed neutrino energy (logarithmic scale) 
for the 107 selected events compared to the MC expectation. 
The neutrino energy is taken to be the sum of the muon mo¬ 
mentum and the energy of any hadronic shower associated 
with the assumed interaction vertex. For the FC and PC 
samples, the muon momentum is determined from range and 
curvature respectively. The solid histogram shows the MC ex¬ 
pectation for the case of no neutrino oscillations, the hatched 
histogram shows the cosmic-ray background and the points 
with error bars show the data. The dashed histogram shows 
the expectation for <-> Vr oscillations with sin^ 2023 = 1-0 
and Amis = 0.0024 eV^ 


for a non-statistical accumulation of events in a particu¬ 
lar region. 

The reconstructed neutrino energy distribution of the 
107 candidate events is shown in FigurelHl The neutrino 
energy is calculated by summing the reconstructed muon 
energy and the hadronic energy of any reconstructed 
shower associated with the start of the muon track. For 
FC events, the muon energy is determined from the track 
range. For PC events the less precise momentum from 
curvature is used. 

The neutrino energy spectrum is sharply peaked to¬ 


wards lower energies and the selected event sample is ex¬ 
pected to have mean neutrino energy of 3.5 GeV (2.0 GeV 
for FG events and 7.0 GeV for the PG events) and mean 
muon energy of 2.4 GeV. For low energy events the abil¬ 
ity to determine the sense of the muon track (up/down) 
is degraded. In MC, 96% of the selected events have the 
correct direction reconstruction. The remaining 4 % of 
events not only have the wrong reconstructed sense, but 
as a consequence, are also assigned the incorrect charge. 
By requiring that |RMSup — RMSdownI > 0.66 ns (see 
Section ns and that the track traverses at least ten 
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planes the fraction of mis-reconstructed events is reduced 
to 0.1%. For the results that follow, the event sam¬ 
ple is divided into two: a ‘Low Resolution’ sample with 
|RMSup — RMSdownI < 0.66 ns and events with ‘Good 
Timing’ for which |RMSup — RMSdownI > 0.66ns. The 
numbers of events in each category are listed in Table^O] 
The 30 events classified as low resolution are mainly short 
events and according to MC have a mean neutrino energy 
of 1.0 GeV. For 85% of the low resolution events, the 
muon is reconstructed with the correct sense (up/down). 
However, in the oscillation analysis that follows, the di¬ 
rection information from the low resolution sample is 
not used, due to the significant fraction of events recon¬ 
structed with the wrong direction sense and the fact that 
for this predominantly low energy sample, the mean an¬ 
gle between the incident neutrino and final state muon is 
large. 


Selection 

Data 

Expected 
no oscillations 

Expected 

Amis = 0.0024 eV^ 

Good Timing 

77 

90 ±9 

68 ±7 

Low Res. 

30 

37 ±4 

28 ±3 

All Events 

107 

127 ± 13 

96± 10 


TABLE III: Classification of events into samples with al¬ 
most unambiguous direction from timing (‘Good Timing’) 
and those where the direction from timing is uncertain (‘Low 
Resolution’). The errors are dominated by the systematic 
uncertainty in the neutrino event rate. The MC expectations 
are given for both no oscillations and sin^ 2023 = 1.0 and 
Amla = 0.0024 eV^ 


Figure cni shows the reconstructed zenith angle dis¬ 
tribution of the 77 candidate events with good timing. 
Of these events, 49 are downward-going (cos0z > 0) 
and 28 are upward-going (cos0z < 0), giving a mea¬ 
sured up-down ratio of 0.571° ;[ 3 (stat.). The statistical 
errors correspond to the 68 % confidence interval calu- 
lated using Poisson statisticsjs^. The expected value 
from Monte Garlo in the absence of neutrino oscillations 
is 0.92±0.03(sj/s.). The expected value is lower than one 
because of the different selection efficiencies for upward- 
and downward-going events and the presence of back¬ 
ground. The upward-going/downward-going double ratio 
is: 

^up‘;dow„/<p%own = 0.62l°l°(stat.) ± 0.02(sys.). 

This is approximately two standard deviations from 
unity, which is the expectation in the absence of neu¬ 
trino oscillations. The systematic error is dominated by 
the experimental uncertainties; the estimated systematic 
uncertainty on t he p redicted up/down neutrino flux ratio 
is less than 1 %|^ due to the relatively high energy of 
the selected neutrino events. 



FIG. 10: The reconstructed cos^zen distribution for the 77 
selected events with good timing compared to the MG expec¬ 
tation. The solid histogram shows the MC expectation for 
the case of no neutrino oscillations, the hatched histogram 
shows the cosmic-ray background and the points with error 
bars show the data. The dashed histogram shows the ex¬ 
pectation for <-> Vr oscillations with sin^ 2023 = 1.0 and 
Amis = 0.0024 eV^ 

A. Oscillation Analysis 

In the two-flavour approximation, which is adequate 
for the level of statistical precision considered here, the 
survival probability, P{v^i —> v^) is given by 

P = 1.0-sin^ 26*23 sin^ ^l-27Am|3[eV^].-^j^^^ , 

where L is the distance traveled by the neutrino and E 
is the neutrino energy. The neutrino path length, L, is 
calculated from the reconstructed zenith angle assuming 
the neutrinos are produced at a height of 20 km in the 
Earth’s atmosphere. Figure ITTI shows the reconstructed 
L/E distribution for the 77 candidate events with good 
timing. 

The reconstructed L/E distribution is used as the ba¬ 
sis for a fit to the hypothesis of Vj- oscillations. 

The resolution on L/E differs greatly event-to-event for 
three main reasons: for PC events the muon momentum 
from curvature may be poorly determined; for low en¬ 
ergy and/or high y events the opening angle between the 
observed muon and the true neutrino direction is large; 
and in the case where the muon direction is close to the 
plane defined by the horizon, relatively small changes 
in angle produce large changes in L/E. To address the 
first issue, PC events with little observable curvature, 
\Q/p\/cr(Q/p) < 1, are not used in the fit to the L/E dis¬ 
tribution. To account for the different L/E resolutions, 
in the oscillation fit the data are binned according to L/E 
resolution. To estimate the event resolution a Bayesian 
approach has been adopted which allows the event-by¬ 
event log (L/E) probability density function (pdf) to be 
determined p^. The rms of the pdf, criogiL/E), gives a 
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FIG. 11: The reconstructed logjQ [L(km)/i5(GeV)] distribu¬ 
tion compared to the expectation. The solid histogram shows 
the MC expectation for the case of no neutrino oscillations, 
the hatched histogram shows the cosmic-ray background and 
the points with error bars show the data. The dashed his¬ 
togram shows the expectation for <-> Vt oscillations with 
sin^ 2623 = 1.0 and Am^g = 0.0024 eV^ 


measure of the log {L/E) resolution of the event. For ex¬ 
ample, Figure II21 shows the data binned in four regions 
of ^\og{L/E)- For MC, the sensitivity of the L/E distri¬ 
bution to neutrino oscillations increases with decreasing 

CTlog(L/£)- 


1. Oscillation Analysis: Fit Procedure 


MINOS 




FIG. 12: The reconstructed logj^Q [L(km)/i5(GeV)] distribu¬ 
tion binned in four regions of log L/E resolution, a. PG events 
with IQ/pI/ujq/p) < 1 are not used. The solid histogram 
shows the MC expectation for the case of no neutrino oscil¬ 
lations, the hatched histogram shows the cosmic-ray back¬ 
ground and the points with error bars show the data. The 
dashed histogram shows the expectation for <-> Vt oscil¬ 
lations with sin^ 2023 = 0.90 and Amlg = 1.3 x 10~®eV^, 
the oscillation parameters corresponding to the best fit to the 
MINOS data. 


The selected events are divided into 10 equal sized 
bins of the estimated uncertainty in reconstructed LjE 
ratio, <J\og(L/E), ranging from 0.1 — 1.1. Events with 
tJiog (L /E) >1.0 are included in the lowest resolution sam¬ 
ple and events with <J\o^(l/e) <0.1 are included in the 
highest resolution. A simultaneous fit is performed to the 
overall normalization (using all selected events), the up- 
down ratio for 77 events with good timing, and separately 
the shapes of the upward and downward LjE distribu¬ 
tions for events with good timing and |Q/p|/ct(q/p) > 1. 
In this way, each event is used only when the physical 
observable being fitted is well-measured. A maximum 
likelihood fit to the data is performed using the negative 
log-likelihood function: 

- In/: = {pi-N\npi)- InP,^) 

k 

- Y,\nft{[LlE],^) - ^ln/,^([L/A],J 

< id 



where N is the total number of observed events and p is 
the total Monte Carlo expectation. The first two terms 


represent the Poisson probability of observing N events 
given the expectation of p. The normalization system- 
atics uncertainties are included as nuisance parameters 
(see below). In the remaining terms, the superscript k 
refers to the bin in criog(L/_E)- The sum 

^(iV^nP„^+Ar,^lnP,^) 

k 

is the “up-down” likelihood. Here and Njl are the 
observed numbers of upward and downward-going events 
with good timing in bin k of resolution; and are 
the Monte Carlo probabilities that an event in resolution 
bin k is upward or downward going {P^ + Pd = !)■ The 
terms 

J2^nf^{[L/E].J and ^ In/|([L/F;].J 

iu id 

are the likelihood functions for the observed LjE distri¬ 
butions of upward- and downward-going events respec¬ 
tively. Here the summations are over the reconstructed 
upward and downward events respectively; /E]i^) 

is the normalized Monte Carlo pdf for the reconstructed 
LjE distribution in the bin of resolution (that of the 
event), evaluated at the measured value of LjE of the 
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event. The MC expectations for /r, P^, P^, f^{L/E) 
and f^{L/E) include contributions from both neutrino 
interactions and cosmic-ray background and depend on 
(Am 23 , sin^ 2 ^ 23 ) and the nuisance parameters represent¬ 
ing the systematic uncertainties. In calculating the ex¬ 
pectations as a function of oscillation parameters, the 
oscillation probabilities are averaged over the distribu¬ 
tion of neutrino production heights obtained from the 
Bartol 3D model. In the fit, systematic effects are in¬ 
cluded using the nuisance parameters, aj, which rep¬ 
resent the deviation of a particular parameter from its 
nominal value. The nuisance parameters contribute to 
the likelihood function through the terms 



where cr^^. is the estimated systematic uncertainty. The 
following systematic effects are included: i) the uncer¬ 
tainty on the expected neutrino event rate is taken to 
be ±10%; ii) a, 3% uncertainty on the muon momentum 
and a 5% uncertainty on the hadronic energy scale; Hi) 
a 3 % uncertainty on the relative efficiency for selecting 
upward- versus downward-going events; iv) to accommo¬ 
date the uncertainty in the shape of the neutrino energy 
spectrum, the spectrum is allowed to scale according to 
1.0+0.ip{E^-2) for E^<2 GeV and 1.0±0.025/3(£;,,-2) 
for Ey > 2 GeV, where (3 is normal distributed (these 
variations cover the differences in the neutrino e nerg y 
spectra obtained from different flux models |^ 1 ^ 1 ^ 1 : 
and v) to allow for uncertainties in the modelling of neu¬ 
trino cross sections, the relative cross section for quasi¬ 
elastic interactions is assumed to have a 20 % uncertainty. 
Since the neutrino flux times cross section is normalized 
to Soudan 2 data, in the fit the systematic error associ¬ 
ated with the QE fraction only affects the shapes of the 
reconstructed LjE distributions. 

With the exception of normalization, the systematic 
uncertainties have little impact on the resulting confi¬ 
dence regions. The above form of the likelihood function 
is chosen to simplify the inclusion of systematic errors in 
the fit: as normalization and up-down ratio are treated 
independently of shape. 

2. Oscillation Analysis: Results 

For each hypothesized value of (Am 23 ,sin^ 26 * 23 ) the 
negative log-likelihood function described above is min¬ 
imized with respect to the nuisance parameters. For 
the data, the minimum likelihood occurs at (Am 23 = 
1.3 X 10“^ eV^, sin^ 2023 = 0.90). The 90% confidence 
limits are obtained from the difference in the log- 
likelihood function —Aln£ = — ln£(Am 23 , sin^ 2 ^ 23 ) ± 
In £0 • Here — In £0 is the value of the negative log- 
likelihood function for the best fit to the data. In the 
limit of Gaussian errors, the 90 % confidence level allowed 
regions of parameter space are defined by — A ln£ < 2.3. 



FIG. 13: The 68 % and 90 % confidence limits on the oscil¬ 
lation parameters obtained using the Feldman and Cousins 
approach. Also shown is the 90 % C.L. limit calculated using 
-Aln£ < 2.3. 


The frequentist approach of Feldman and Gousinsls^ is 
used to determine the value of — Aln£ which corresponds 
to a particular confidence level. For each point in param¬ 
eter space, (Amgg„,sin^ 26gen), 1000 Monte Carlo exper¬ 
iments are generated. In each experiment a value for 
each systematic bias is drawn from a Gaussian distri¬ 
bution with standard deviation equal to the estimated 
systematic uncertainty. The generated data sample is 
fitted in the same manner as the data and the value of 
—Aln£(Am^g„, sin^ 29gen) is determined. The generated 
point is included in the 90 % confidence region if less than 
90% of the experiments yield a smaller value than ob¬ 
tained in the data. For the MINOS data the 68 % and 
90% confidence limits obtained using the Feldman and 
Cousins approach are shown in FigureCSl The 90% C.L. 
contour is close to that obtained using — Aln£ = 2.3. 
With the current statistics, the MINOS atmospheric neu¬ 
trino data are consistent with a wide range of oscillation 
parameters including the most recent results from Super- 
Kamiokande0 and K2KpiT|. The data disfavor the null 
oscillation hypothesis at the 98 % confidence level. 

For completeness. Figure ITU shows the likelihood as a 
function of Am 23 for sin^ 2023 = l-O- The rises at large 
and small values of Am |3 are mainly due to the nor¬ 
malization and up-down ratio. The structure within this 
broad minimum arises from the fit to the shape of the 
LjE distribution. The quality of the fit is good. As a 
measure of the fit quality, 10000 simulated experiments 
were generated with (Am 23 = 0.00 24eV^, sin^ 2023 = 
I.O) and the minimum value of —ln£o determined; in 
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FIG. 14: The A/! curve as a function of Am^s for the case 
of maximal mixing. The dotted curved shows the expected 
sensitivity which is the average likelihood curve obtained from 
10000 Monte Carlo experiments. 

84% of these experiments the minimum value of —In/lo 
exceeded that obtained from the fit to the data. FigurellH 
also shows the expected sensitivity. 

B. Charge Ratio 

The selected contained events with unambiguous tim¬ 
ing information are divided into neutrino and anti¬ 
neutrino interactions on the basis of the reconstructed 
muon charge obtained from the curvature of the recon¬ 
structed muon track. Only events with unambiguous di¬ 
rection from timing are considered; events with an in¬ 
correct direction will be reconstructed with the wrong 
zenith angle and will have their charge inverted. Fig¬ 
ure CHI shows the distribution of {Qlv)l for the 77 
events with well determined direction from timing com¬ 
pared to the MC expectation. Events are classified as: 

for < -2; for (Q/p)/cr(Q/p) > +2; 

or events which are ambiguous, \{Q/p)/( q/ p)\ < 2. The 
selected numbers of events in each charge category are 
compared to the MC expectations in Table CYI 

Of the events where it is possible to cleanly tag the 
charge of the muon, 18 are identified as Vp, candidates 
and 34 as candidates, yielding a 77^ to ratio of: 

= 0.53t°:?i(stat.) ± 0.03(sys.). 

The systematic uncertainty is the experimental uncer¬ 
tainty associated with charge identification. The uncer¬ 
tainty was estimated by shifting and smearing the Monte 
Carlo reconstructed values of {Q/p)l(J{Q/p) whilst main¬ 
taining reasonable agreement between data and Monte 
Carlo for the stopping muon data shown in FigureOl For 
the purposes of studying possible biases in the charge 
reconstruction, 37 % of the data were recorded with the 
coil current reversed. Consistent values for the 



FIG. 15: The reconstructed distribution of (Qlp)l(y(Q/p-j, the 
ratio of the charge divided by momentum obtained from the 
track curvature divided by its error. The solid histogram in¬ 
dicates the Monte Carlo expectation assuming no oscillations, 
the hatched histogram shows the cosmic-ray background and 
the points with error bars show the data. The dashed his¬ 
togram shows the expectation for Vp ^ Vr oscillations with 
sin^ 26123 = 1.0 and Am^g = 0.0024 eVT 


Selection 

Data 

Expected 
no oscillations 

Expected 

Amis = 0.0024 eV^ 

Low Res. 

30 

37 ±4 

28 ±3 

Ambig. Vp,IVp 

25 

26 ±3 

20 ±2 

Vp 

34 

42 ±4 

31 ±3 

Vp. 

18 

23 ±2 

17 ±2 


TABLE IV: Event classihcation according to timing and track 
curvature. The four categories are events with ambiguous 
direction from timing (‘Low Resolution’), events with good 
timing information but ambiguous charge assignment (‘Am¬ 
biguous Vp./Vfi'), Vp and Vp. The errors are dominated by the 
systematic uncertainty in the neutrino flux x cross section. 
The MC expectations are given for both no oscillations and 
sin^ 26123 = 1.0 and Amig = 0.0024 eVT 

ratio are found in the normal and reversed current data 
samples; 22(12) events are identified as and 12(6) are 
identified as in the normal(reversed) field data sam¬ 
ples. From Monte Carlo, the expected ratio of identified 
Vp to Vp events is 0.550, where it is assumed that both 
neutrinos and anti-neutrinos oscillate with the same pa¬ 
rameters. This expected ratio is almost independent of 
the values of the oscillation parameters provided they 
are the same for neutrinos and anti-neutrinos. The ratio 
of Vp to Vp events in the data compared to the Monte 
Carlo expectation (Bartol 3D and NEUGEN3) assuming 
the same oscillation parameters for neutrinos and anti¬ 
neutrinos is: 

= 0.96™(stat.)±0.15(sys.). 

The statistical errors correspond to the 68 % confidence 
interval calulated using Poisson statistics|3^. The sys¬ 
tematic error includes the experimental uncertainty as- 
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sociated with the muon charge identification (0.06), the 
uncertainty in the relative 17^ to flux (0.04) and the rel¬ 
ative uncertainty in the to cross section (0.13). The 
systematic errors on the relative to fluxes and cross 
sections were estimated taking into account the energy 
spectra of the charge-tagged neutrino and anti-neutrino 
events. From Monte Carlo, the sample of events where 
the charge of muon is cleanly identified has a mean neu¬ 
trino energy of 3.7 GeV, with 95 % of the expected events 
having neutrino energies between 1 GeV and 10 GeV. In 
Monte Garlo, 40 % of these events arise from quasi-elastic 
interactions, 30 % arise from resonance production, and 
30 % from deep inelastic scattering. In this energy range 
the uncertainty on the ratio of atmospheric to flux 
was estimated to be 8%[^ for the Bartol ID model. 
Recent studies based on the Bartol 3D model give an 
estimated uncertainty of 4%[^. Since there is limited 
data on anti-neutrino cross sections in the energy range 
1 — 5 GeV [ 33 ,m, 113, the variation in the predicted 
event rate was studied through conservative changes to 
the neutrino cross section model. The variations consid¬ 
ered include parameters affecting the free nucleon cross 
sections such as axial vector masses, choice of PDF set, 
and model for the resonance region. Similarly the effect 
of changes to the nuclear physics model which affect the 
rate via Pauli blocking of quasi-elastic and nuclear shad¬ 
owing of DIS events was studied. The quadrature sum 
of these changes is 13.5%. The largest contribution to 
this uncertainty comes from the treatment of resonance 
production and the resonance/DIS transition region. A 
12 % difference is found in comparing a model which ex¬ 
plicitly includes resonance producti on l4Il| versus one 
which uses a QCD-based approach The size of this 
difference ultimately reflects the uncertainty in the ex¬ 
perimental data to which these models are tuned. 

CPT violating models which attempt to explain the 
LSND data SMgest a large value of Am^g for anti¬ 
neutrinos [ 13 . In principle, the MINOS data will be 
able to address this possibility by measuring the oscil¬ 
lation parameters for the selected anti-neutrino sample. 
Figure M shows up/down distribution of the 18 and 
34 events compared to the expectation for: i) no oscil¬ 
lations; ii) the case where both and oscillate with 
Am^g = 0.0024 eV^ (maximal mixing) and; iii) the case 
where oscillate with Am^g = 0.0024 eV^ (maximal 
mixing) and F^ oscillate with Am^g = 1.0 eV^ (maximal 
mixing). The data are consistent with the same oscilla¬ 
tion parameters for neutrinos and anti-neutrinos. How¬ 
ever, with the current statistics the possibility of a large 
value of Am 2 g for anti-neutrinos cannot be excluded. 


VII. SUMMARY 

The MINOS far detector has been taking data since 
the beginning of August 2003 at a depth of 2070 me¬ 
ters water-equivalent in the Soudan mine, Minnesota. 
This paper presents the first MINOS observations of 


and F^ charged current atmospheric neutrino interactions 
based on an exposure of 418 days. A total of 107 candi¬ 
date contained-vertex neutrino interactions is observed, 
consistent with both the expectation of 127 ± 13 for no 
neutrino oscillations and 96 ± 10 for Amfg = 0.0024 eV^ 
and sin^ 2023 = 1-0. The expected numbers of events in¬ 
clude the estimated background from cosmic-ray muons, 
4.4 ± 0.5, obtained from data. The errors on the expec¬ 
tation are dominated by 10 % uncertainty on the neu¬ 
trino event rate which was obtained using results from 
the Soudan 2 collaboration. Of the events for which the 
direction can be cleanly identified, the ratio of upward 
to downward-going events in the data is compared to the 
Monte Carlo expectation in the absence of neutrino os¬ 
cillations giving: 

^np/dow„/<p/down = 0± 0.02{sys.). 

An extended maximum likelihood fit to the ob¬ 
served log L/E distribution yields a best fit value 
of (Amfg = 1.3 X 10“^ eV^, sin^ 2023 = 0.90) and 
90% confidence limits of (7 x 10“^ eV^ < Amfg < 
5 X 10“^ eV^, sin^ 2023 > 0.2). The consistency of the 
data with the null hypothesis of no neutrino oscillations 
is investigated; the data exclude the null hypothesis at 
the 98 % confidence level. 

The curvature of the observed muons in the 1.3 T MI¬ 
NOS magnetic field is used to separate and F^ inter¬ 
actions. Of the selected events for which it is possible to 
cleanly determine the charge of the muon, 18 are identi¬ 
fied as F^ candidates and 34 as Vfj, candidates, giving an 
observed F^ to ratio of 0.53~!iQ'il{stat.) ± 0.03(sj/s.). 
The fraction of F^ events in the data is compared to the 
Monte Carlo expectation assuming neutrinos and anti¬ 
neutrinos oscillate in same manner giving: 

= 0.96™(staf.)±0.15(sj/s.). 

Although the statistics are limited, this is the first di¬ 
rect observation of atmospheric neutrino interactions sep¬ 
arately for z/^ and F^. The data are consistent with neu¬ 
trinos and anti-neutrinos oscillating with the same pa¬ 
rameters, although CPT violating scenarios with large 
values of Am|g for anti-neutrinos are not excluded with 
the current data. 
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FIG. 16: The charge-separated up/down distributions of the events which are cleanly identified as or V^. The data (r'^ on 
the left and 77^ on the right) are compared to the expectation for i) no oscillations; ii) the case where both and 77^^ oscillate 
with Am 23 = 0.00 24 eV^ and; iii) the case where oscillate with Am |3 = 0.0024 eV^ and 77^ oscillate with Am 23 = 1.0 eV^. 
For the oscillated predictions maximal mixing is assumed for and 77^. 
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